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ABSTRACT 

A Reduced Xavier-Stokes solution technique has been sucessfully used to design vortex 
generator installations for the purpose of minimizing engine face distortion by restructuring the 
development of secondary flow that is induced in typical three-dimensional curved inlet ducts. 
The results indicate that there exists an optimum axial location for this installation of co-rotating 
vortex generators, and within this configuration, there exists a maximum spacing between gen- 
erator blades above which the engine face distortion increases rapidly. Installed vortex generator 
performance, as measured by engine face circumferental distortion descriptors, is sensitive to 
Reynolds number and thereby the generator scale, i.e. the ratio of generator blade height to local 
boundary layer thickness. Installations of co-rotating vortex generators work well in terms of 
minimizing engine face distortion within a limited range of generator scales. Hence, the design 
of vortex generator installations is a point design, and all other conditions are off-design. In 
general, the loss levels associated with a properly designed vortex generator installation are very 
small, thus they represent a very good method to manage engine face distortion. This study also 
showed that the vortex strength, generator scale, and secondary flow field structure have a com- 
plicated and interrelated influence over engine face distortion, over and above the influence of the 
initial arrangement of generators. 


INTRODUCTION 

Modern tactical aircraft are required to be maneuverable at subsonic, transonic, and 
supersonic speeds, without giving up good cruise performance. Consequently, proper integration 
of the engine inlet with the airframe is of paramount importance. Regarding the enhancement 
of inlet performance and operation, design for optimum airframe-inlet integration has the fol- 
lowing goals: (1) to minimize approach flow angularity with respect to the inlet cowl lip, (2) to 
deliver uniform, high pressure recovery flow to the inlet face, (3) to prevent or minimize vortex, 
wake, and boundary layer ingestion by the inlet throughout the flight envelope (4) to reduce 
FOD/hot gas ingestion by the inlet, and finally (5) to minimize the potential for flow field inter- 
ference from weapon carriage, 'firing, landing gear deployment, tanks, pods, or other hardware. 
The combination of inlet design and airframe integration must not only provide high pressure 
recovery to maintain the desired thrust levels, but also generate low flow distortion consistent 
with stable engine operation. 

Engine face flow distortion is one of the most troublesome and least understood problems 
for designers of modern inlet engine systems (Refs. 1 and 2). One issue is that there are numerous 
sources of flow field distortion that are ingested by the inlet or generated within the inlet duct 
itself. Among these sources are (1) flow separation at the cowl lip dunng maneuvenng flight, (2) 
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flow' separation on the compression surfaces due to shock-wave boundary layer interactions, (3) 
spillage of the fuselage boundary layer into the inlet duct, (4) ingestion of aircraft vortices and 
wakes emanating from upstream disturbances, and (5) secondary 7 flow and possibly flow sepa- 
ration within the inlet duct itself. Most aircraft have experienced one or more of these types of 
problems during development, particularly at high Mach numbers and/or extreme maneuver 
conditions, such that flow distortion at the engine face exceeded allowable surge limits. Such 
compatibility problems were encountered in the early versions of the B70, the F-lll, the F-14, 
the MIG-25, the Tornado and the Airbus A300 to name a few examples. 

One of the most commonly used methods to control local boundary layer separation 
within inlet ducts entails the placement of vortex generators upstream of the problem area. 
Vortex generators in use today are small wing sections, mounted on aircraft wing surfaces or the 
inside surfaces of inlet ducts, inclined at an angle to the oncoming flow to generate a shed vortex. 
The generators are usually sized to local boundary layer height to allow for the best interaction 
between the shed vortex and boundary layer itself, and are usually placed in groups of two or 
more upstream of the problem area. The principle of boundary layer control by vortex generators 
relies on the induced mixing between the external or core stream and the boundary' layer region. 
This mixing is promoted by vorticies trailing longitudinally over the duct surface adjacent to the 
edge of the boundary’ layer. Fluid particles with high momentum in the streamwise direction are 
swepted along helical paths towards the duct surface to mix with and, to some extent, replace the 
low momentum boundary layer flow. This is a continuous process that provides a source of re- 
energization to counter the natural boundary' layer growth caused by friction, adverse pressure 
gradients, and low energy secondary flow accumulation. There are two basic configurations of 
vortex generators. In one configuration, all the vortex generators are inclined at the same angle 
with respect to the oncoming flow direction. These are called co-rotating configurations because 
the shed vorticies rotate in the same direction. In the other configuration, the vortex generators 
are grouped in pairs inclined in the opposite direction to the flow, such that pairs of counter- 
rotating shed vorticies are generated. 

Co-rotating vortex generators are very competitive in reducing flow’ separation if the gen- 
erators are properily selected and located. The main advantage of co-rotating type vortex gen- 
erators are their downstream effectiveness, resulting in more effective usage of the vortex energy 
w’ithin the affected boundary layer. According to design "wisdom", this type of vortex generator 
has a few’ special advantages when used w’ithin S-duct inlet configurations, namely: (1) the in- 
duced vortices will remain close to the W’all; consequently a "cleaner" core flow’ w’ill result, and (2) 
the induced vorticies will counteract the natural and often strong secondary flow’s w’hich can de- 
velop. 

Counter-rotating, equal strength vortex generators have been used in a number of aircraft 
inlet ducts, the F/A-18 and the center inlet duct on the production 727 aircrafts, to name a few’ 
examples. This type of vortex generator is very’ effective in reducing flow’ separation if the vortex 
generators are placed slightly upstream of the region of separation. How'ever, according to vortex 
generator design "wisdom", the disadvantages of this type of generators, as compared to co- 
rotating generators, are that the induced vorticies tend to lift off the duct surface, thus reducing 
their effectiveness, causing higher loss in inlet recover)’ and larger total pressure distortion at the 
compressor face. 

The performance of vane-type vortex generators was evaluated by Taylor (Ref. 3), for 
diffusers and airfoils at low speed, and by Valentine and Carrol (Refs. 4 and 5), for airfoils and 
wings at high speeds. This work provides trends in effectiveness for certain vortex generator de- 
sign variables such as angle-of-attack, height, distance ahead of separation, etc. Attention was 
focused on the detailed changes that were produced in the boundary layer as a result of placement 
of vortex generators in the flow’. Percy and Stuart (Ref. 6), extended the study of the effects of 
various design parameters and concluded that the strength and disposition of the individual in- 
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duced vortices was more important than the details of the boundary layer upstream of the im- 
posed pressure gradient. 

In general, the design strategy adopted in these studies was basically a boundary' layer 
approach, whereby the vortex generators provided a continuous re-energization of the boundary 
layer to prevent flow separation. Thus, improved engine face distortion levels were achieved by 
preventing local flow separation within the inlet duct. However, vortex generators within inlet 
ducts did not always perform as expected, primarily because the flow environment in which they 
were positioned was not properly taken into consideration. The "design rules" for the geometry 
and placement of vortex generators were based on two-dimensional boundary layer concepts, and 
they did not operate well in regions with high secondary flow. 

It was not until the confirmation test by Kaldschmidt, Syltebo, and Ting (Ref. 7), for the 
refanned JT8D engine on the 727 center duct inlet that an attempt was made to use vortex gen- 
erators to restructure the development of secondary flow in order to improve the engine face 
distortion level. Thus, a very important shift in strategy on the use of vortex generators had oc- 
curred. The perspective had moved from a local two-dimensional boundary layer approach to 
eliminate local flow separation, to a global three-dimensional vortex-secondary flow interaction 
concept, where the design goal was now to control the development of three-dimensional sec- 
ondary flow itself. 

In order to accomplish this new objective for internal flow control, the design strategy 
must shift from an experimental to an analysis based methodology. A definitive approach for 
computation of three-dimensional flow in curved inlet passages with vortex generator control 
would be the numerical solution of the full Navier-Stokes (FNS) equations. Although feasible, 
in three-dimensions this approach requires considerable computational effort, especially in view 
of the necessity of using a multi-zone strategy to resolve the very small length scales associated 
with the vortex generator region along with the much larger length scale phenomena associated 
with the inlet duct itself. A more cost effective method for predicting three-dimensional turbulent 
subsonic vortex flows in curved passages would be solutions to the reduced form of the Navier- 
Stokes (RNS) equations. The objective of this approach is to introduce approximations which 
adequately represent the essential physical processes of. interest and yet lead to governing 
equations which can be solved much more economically than the full Navier-Stokes equations. 

This paper represents one in a series of studies on the design issues associated with inlet- 
engine compatibility problems. These studies center on the development of CFD tools and 
techniques which look promising within an analysis-design environment, and the application of 
these new analysis approaches to understanding and controlling inlet-engine distortion. The first 
paper in this series by Anderson (Ref. 8), deals with the aerodynamic characteristics of vortex 
interaction within the F/ A- 18 inlet duct, where the vortex interaction arises as a result of a vortex 
ingestion. Later studies will envolve the effect of Vortex ingestion on the engine face flow field 
itself. 


The overall goal of this effort is to advance the understanding and control of engine face 
distortion, and in particular, to analyze the basic interactions that can influence this important 
design problem. Specifically, the current paper achieves two goals, namely: (1) the development 
and validation of subsonic RNS computational techniques for analysis of general geometry inlet 
ducts with "real" vortex generator model capabilities, and (2) the formulation of a design strategy 
using CFD for the use of vortex generators to control the development of secondary flow within 
inlet ducts to minimize the engine face distortion. 


3 



THEORETICAL BACKGROUND 

The reduced Xavier-Stokes (RXS) equations, originally termed parabolized Navier-Stokes 
(PX'S) equations, and more recently semi-elliptic, or partially parabolic, are used here as an 
initial-value space marching method for the evaluation of subsonic compressible flow with strong 
interactions and/or separation arising from internal vortex flows. Techniques that use space 
marching with an approximate form of the RN T S equations, namely initial-value methods and 
those that require three-dimensional global iterations, have been used for a number of years to 
predict flows in curved ducts and turbomachinery blade cascades. Unfortunately, this terminol- 
ogy does not identify the relevant mathematical approximations nor does it distinguish these ap- 
proximations from the properties of the solution algorithm and the differential or difference 
equations. In other words, different methods within the same '’category" will in some instances 
give significantly different results. 

Partitioned Geometry and Mesh Generation 

Three-dimensional viscous subsonic flows in complex inlet duct geometries are investigated 
by a numerical procedure which allows solution by spatial forward marching integration, utilizing 
flow approximations from the velocity-decomposition approach of Briley and McDonald (Refs. 
9 and 10). The goal of this approach is to achieve a level of approximations that will yield ac- 
curate flow predictions, while reducing the labor of solution below that of the full X’avier-Stokes 
equations. The governing equations for this approach have been given previously for orthogonal 
coordinates, and the approach has been applied successfully to problems whose geometries can 
be fitted conveniently with orthogonal coordinate systems (Ref. 11). However, geometries en- 
countered in typical subsonic inlet ducts cannot be treated easily using orthogonal coordinates, 
and this lead to an extension of this approach by Levy, Briley, and McDonald (Ref. 12), to treat 
ducted geometries with nonorthogonal coordinates. The nonorthogonal capability has been val- 
idated over a wide range of inlet flow conditions by Towne (Refs. 13 and 14), and Anderson (Ref. 
15). 


Although the analysis itself was general, the class of ducted geometries that could be ana- 
lyzed was represented by superelliptic cross-sections normal to a reference line space curve having 
continuous second derivatives. The description of the superelliptic cross-sections was specified 
by polynomials defined in terms of a marching parameter, t. In generalizing the geometry for- 
mulation Anderson (Ref. 8), extended the analysis to cover ducted geometries defined by an ex- 
ternally generated gridfile. This version of the 3D RNS computer code is called RNS3D. The 
geometry description within the gridfile is a "ducted" geometry which has a variable cross- 
sectional area and shape and a centerline which is curved and possibly twisted. In addition, the 
duct described by the gridfile is considered to have a defined centerline with continouus second 
derivatives. The surface geometry is described in terms of cross-sectional planes which lie per- 
pendicular to the duct centerline, and thus represent the flow area at each streamwise station. 
Since the inlet duct geometry definition has been reduced to a cross section specification which 
is placed perpendicular to a centerline space curve, then a number of grid and geometry’ pre- 
processing functions may be performed using RX T S3D. These pre-processing functions include: 
(1) reclustering the existing gridfile mesh points distribution for more accurate solutions in re- 
gions of high shear, (2) redefining the centerline space curve to satisfy design constraints, and (3) 
altering the cross-sectional shape of the inlet duct to reflect specified design iterations. The ap- 
proaclf taken by Anderson (Ref. 3), is to develop a geometry pre-processor to augment the ex- 
isting geometry and grid generation programs for internal inlet duct configurations, i.e. to 
partition the "work" of mesh generation between the grid generator and flow solver. 
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Vortex Generator Model 


The model for the vortex generators within the RNS analysis, described by Kunik (Ref. 
16), takes advantage of the stream function-vorticity formulation of the governing equations. 
The shed vortex is modeled by introducing a source term into the vorticity equation that is a 
function of the geometric characteristics of the generators themselves. This source term is in- 
troduced at every' point in the cross-plane in the form of the following expression 

t p = rv" (c,f2) (i) 

where T, is the vortex strength at any point in the crossplane, To is the vortex strength at the tip 
of the generator, r is the distance between the field point and at the tip of the generator, and C\ 
is a constant which controls the decay of the shed vortex strength in the crossplane. This vortex 
model resembles the one proposed by Squire, (Ref. 17), except that it neglects the variation of 
viscosity in the cross-plane. Downstream of the generator region, the analysis uses the turbulence 
model of McDonald and Camarata (Ref. 18), without adjustments for the influence of the gen- 
erators themselves. 


The FLARE Approximation 

The analysis as presented here is applicable only when the primary velocity is not negative. 
Since "small" regions of reverse flow can arise in curved inlet ducts, the numerical method is lo- 
cally' modified to permit forward marching when the flow contains small regions of reverse flow. 
The technique used follows Reyhner and Flugge-Lotz approximation, (Ref. 19), by adding small 
artificial convection at grid points where the primary flow is reversed. This is known as the 
FLARE approximation, after the authors of Ref. 19. For thin regions of reverse flow, although 
the area of flow separation can be very large compared to the passage itself, the technique permits 
the analysis to proceed downstream beyond reattachment, confining the FLARE approximation 
to the separated region. 


Steady State Engine Face Distortion Descriptors 


It is impractical to measure anything at the engine face when the engine is installed and 
operating, consequently, the engine and inlet designers agreed upon an Aerodynamic Interface 
Plane (AIP) which is forward of the compressor face but sufficiently close to the engine face to 
have a similar flow field. Current L.S. practice uses forty or forty-eight transducer probes ar- 
ranged in eight rakes with five or six rings. The radius of each ring is set such that all probes are 
at the centroid of equal areas. All distorton descriptors, whether they quantify steady state or 
transient distortion conditions, are always calculated relative to the standard rake located at the 
AIP. 

The most widespread quantitative distortion descriptor available in the literature, because 
of its use in the earliest measurements on inlet ducts in the late 1950's, is simply: 


Dt = 


Pt —Pi 

1 *max 1 *min 

Pi 

1 l a\e 


( 2 ) 


where Pi miX is the maximum rake total pressure, Pt mn is the minimum rake total pressure, and 
Pi mi is the area weighted average rake total pressure. In experimental data reduction, it is as- 
sumed that the both the static pressure and temperature are constant and steady across the 
aerodynamic interface plane (AIP); thus both the velocity and Mach number can be considered 
functions only of total pressuc and the distribution of this quantity is the only measurement that 
needs to be made. This parameter is always useful to determine for comparison purposes and to 
describe the 'general health' of inlet ducts irrespective of the type of powerplant that may be used. 
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The effect of circumferential distortion on compressor surge margin is essentially to drop 
the maximum pressure ratio of a constant corrected speed line. One of the simplest quoted de- 
scriptor for circumferential distortion is from Rolls Royce and is defined as 

Pt -Pt 

I i ave l < mm 

tfave 

where Pt„, and q m . are the average total and dynamic pressure at the engine face or aerodynamic 
interface plane and Pt mm is the mimimum total pressure in any pi-section of extent 6. Significant 
6 values can vary with engine design and commonly are 60°, 90°, and 120°. For bypass engines, 
a circumferential distortion descriptor DC,-. C c is often used, where GG indicates that the index 
is taken over the area of the gas generator. 




More advanced distortion descriptors, introduced in the late 1960's and 1970's, take into 
account the Dt distortion of each ring of total pressure measurements. Thus, the radial distortion 
Dt, is defined as 


Dl r 


Pi - Pt 

1 *max 1 l ave 

Pt 

* * may 



(4) 


where Pt m , is the average total pressure for a given ring radius and Pt ma 
total pressure. The circumferential distortion Dt e is defined as: 


Dig — 


Pt -Pt 

1 *gvg * *mm 

Pt 



is the maximum local ring 


(5) 


where Pt mm is the lowest average total pressure in any 6 segment, usually 60° or 180° of arc for 
a given ring radius having an average ring total pressure Pt„,. The ring distortion descriptors 
Dt, and Dt e are both functions of engine face radius, and will be expressed in terms of normalized' 
engine face radius, i.e., 


( r cowl- r dome ) ^ 

where r is the engine face field point radius, r dom , is the radius of the engine dome, and is the 
engine cowl radius. 


RESULTS AND DISCUSSIONS 
Comparison and Validation with Experimental Data 

To demonstrate the accuracy of the numerical results obtained with RNS3D for internal 
duct flows typical of high angle-of-attack conditions, a series of numerical simulations -were car- 
ried out using the University of Tennessee diffusing S-duct. In this experimental investigation 
sponsored by NASA Lewis Research Center Vakili, Wu, Liver, and Bhat (Ref. 20), obtained a 
series of measurements in a 30-30 degree diffusing S-duct of area ratio 1.5 with and without 
vortex generators (Fig. 1). The 30-30 degree circular cross-section S-duct, shown in Fig. 1, was 
made from two symmetric sections. The inlet duct diameter D, was 16.51 cm. and duct cenerline 
had a mean radius of curvature R of 82.55 cm. A straight pipe section of length 4.7 5D, was in- 
stalled upstream of the curved section to allow for the development of the turbulent boundary 
layer to the desired thickness. Another pipe of length 9.0 D, was installed downstream of the S- 
duct. 


All measurements were made at a nominal inlet Mach number of 0.60 at the reference 
measurement station in the straight section at XI D, = — 1.54. The flow parameters at this station 
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were used as reference conditions for non-dimensionizing the experimental data. The exper- 
imental survey stations correspond to X/A = 0.0, 1.29, 2.49, and 5.2. At each survey station, a 
five-port cone probe was traversed radially at ten azimuthal angles, approximately 20 degrees 
apart, on both sides of the symmetry plane. At least seventy points were measured at each 
traverse. 

A polar grid topology (Fig. 2) was chosen for the University of Tennessee diffusing S-duct, 
consisting of 49 radial, 49 circumferential, and 101 streamwise nodal points in the half-plane. 
The internal grid was constructed such that the transverse computational plane was perpendicular 
to the duct centerline. Grid clustering was used both in the radial and circumferential directions 
to redistribute the nodal points along these coordinate lines to resolve the high shear region near 
the wall and the separation region in the second bend. The flow in the duct was turbulent, with 
an entrance Mach number of 0.6, Reynolds number based on hydraulic diameter of 1.76^10*, and 
a shear layer thickness of <5 '£>, = 0.05. These initial conditions were applied at approximately an 
axial station 1.54 inlet diameters (D,) upstream of the duct entrance. 

Fig. 3 shows a comparison between the experimental and computed total pressure coeffi- 
cient contours at X; D,= 5.2. In both the experiment and analysis, the flow in the S-duct sepa- 
rated and reattached in the second bend upstream of XID,= 5.2. This flow separation was caused 
by both adverse streamwise pressure gradient and the effect of pressure-driven secondary flow 
resulting from duct curvature. Experimental measurements and computational results from the 
3D RXS code include the phenomena of separation and reattachment upstream of X/Z>,= 5.2 and 
show excellent agreement for a simple mixing length turbulence model. 

As separation was encountered in the second bend of the S-duct, three pairs of vortex 
generator devices were installed in the duct at X/Z),=0.09, and at circumferential angles of -38.0, 
0.0, and + 38.0 degrees. The vortex generator pairs had geometric incidence angles of + 16.0 and 
-16.0 degrees relative to the duct centerline. Fig. 4 shows the comparison between the exper- 
imental and computed total pressure coefficient contours at X/D, = 5.2. Comparison of contour 
levels between the separated case (Fig. 3), and the vortex generator case (Fig. 4), shows that the 
vortex generators successfully mixed the high energy core flow with the low energy flow in the 
wall region to suppress separation. In general, the computed interaction between the induced 
vortex generator flow and the pressure driven secondary flow was physically realistic and the 
agreement between experiment and analysis is considered very good, although more improve- 
ments on the generator model must be made. 

Figs. 5 and 6 show additional flow characteristics obtained with the 3D RNS analysis with 
the vortex generator modeling. The secondary flow structure from the vortex generator model 
just downstream of the generator region, i.e. at XI D, = 0.18 is shown in Fig. 5, and clearly reveals 
the three pairs of vortices that arise from the three pairs of counter-rotating generators. The 
limiting streamline signature shown in Fig. 6 indicates that the generator configuration tested 
eliminated the flow separation encountered in the second bend, and reveals the familiar 
topographical pattern through the generator region itself. 

Flow Separation and Vortex Liftoff 

The three-dimensional separation encountered in the University of Tennesse diffusing S- 
duct was very large in area and thin in extent (Figs. 7 and 8), and this separation did not altered 
the pressure distribution in a substantial manner. Secondary flow resulting from duct curvature 
caused an accumulation of boundary layer near the innerwall of the first 30 degree bend. The 
thick boundary layer thus established was especially susceptible to flow separation because of the 
adverse streamwise pressure gradients induced by the reverse curvature section of the second 30 
degree bend. A comparison between the computed oil flow patterns (represented by the limiting 
streamline topology in Fig. 7) and the experimental oil flow patterns presented in Fig. 8 shows 
excellent correspondence. Of exceptional importance is the fact that the space marching RXS 
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analysis method using FLARE approximations captures the reverse flow region of this sepa- 
ration. 


To meet the required confidence level of code validation, it is also essential that the 3D 
RNS marching analysis be able to capture the known topological structure of the limiting 
streamlines in the vicinity of a three-dimensional separation. A very striking and significant fea- 
ture captured by the analysis (Fig. 9), and seen in the oil flow pattern (Fig. 8), is the convergence 
of the limiting streamlines as an indication of three-dimensional separation taking place in this 
duct. Another important and striking feature is the symmetric pair of spiral nodes and pair of 
saddle points that were clearly captured by the 3D RNS analysis. The topological patterns, as 
shown in the analysis Fig. 9, and the photograph of the surface oil flow pattern Fig. 8, also reveal 
the remarkable characteristic that the limiting streamlines forming the spiral node enter only from 
downstream of the nodal point. • The very familiar topological pattern shown in Figs. 7 and 8 is 
known to describe the important stage in the development of the pair of counter-rotating vortices 
that form in the first 30 degrees of turning resulting in vortex liftoff in the second bend. 

The Design Problem Using Vortex Flow Control 

Unlike external aerodynamics where the two-dimensional flow assumption is valid over 
extended regions of the flow field, the internal flow problem is plagued with three-dimensional 
flow effects. When viscous effects or other sources of vorticity are present, three-dimensional 
flows differ fundamentally from their two-dimensional counterparts in that large secondary flows 
are generated by a deflection of the primary flow and/or other mechanisms. Secondary flow 
theory (reviewed by Horlock and Lakshminaryana Ref. 21, and Lakshminarayana and Horlock 
Ref. 22) affords considerable insight into the generation of secondary flow and establishes that 
large secondary flows can be generated by small deflections of vorticity or shear. The large sec- 
ondary flows thus generated often exert an appreciable influence on the primary flow, and thus 
aerodynamic performance, viscous losses, and engine face distortion can be significantly affected. 
Therefore, the design problem is to control the three-dimensional secondary flows that are gen- 
erated within typical inlet ducts through an arrangement of vortex generators for the purpose of 
minimizing engine face distortion, and the effectiveness of the the design is judged by standard 
engine distortion descriptors. 

Inlet Duct and Vortex Flow Control Design Variables 

The 727, JT8D-100 center inlet duct geometry and computational grid used in this study 
is shown in Fig. 10, and the single block polar cross-sectional grid topology, with 99 radial and 
49 circumferential grid points in the half-plane is presented in Fig. 10 at two axial stations, i.e. 
at XI R, = 0.0 and X/R, = 5.0. The computation were made at an inlet enterence Mach number 
of 0.6, Reynolds number of 12-OxlO 6 based on hydraulic inlet diameter (D,), and a shear layer 
thickness 5/ D, = 0.005. 

The geometry of the co-rotating vortex generators used in this study along with the no- 
menclature used in positioning the individual blades are presented in Figs. (11) and (12). The 
important geometric design parameters include: (1) the vortex generator blade height ( h/R ,), (2) 
the blade chord length (c/R,), and (3) the vane angle of attack (/?*). For all the calculations within 
this study, the vortex generator blade height ( hjR ,) was set a 0.075, the ratio of generator height 
to chord length (h/c) was fixed at 0.5, and the vane angle of attack (/?„) was set at 16.0°. Instead 
of the usual spacing parameter ( d/R ,), i.e. the distance between adjacent blades, the positioning 
of the vortex generator blades was described in terms of a spacing angle (a.j) and a sector angle 
over which the blades were positioned (#,). For this study, the relationship between blade spacing 
angle (y., s ) and sector angle ( d s ) is given by 
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where (nj is the number of vortex generator blades. Equ. (7) was also used to position the in- 
dividual generator blades around the inside periphery of the inlet duct at a given axial sector lo- 
cation ( XJR ,). The angle (d s ) was measured counter-clockwise relative to an azimuthal angle of 
180° with respect to the vertical axis of the duct. 

Shown in Fig. 13 are the axial locations of the vortex generator sector regions (XJR,) 
covered in this study. These sector regions were located between XJRi = 2.0 and XJRi = 7.0, and 
cover a sector angle (0,) up to 157.5° as measured counter-clockwise relative to an azimuthal 
angle of 180° with respect to the vertical axis of the duct. 

Tables I, II, and III present the range of vortex generator configurations, or test cases, 
considered for this study, and these have been arranged in five different groupings. In the first 
grouping, i.e. Configs. 1 through 4, the effect of vortex generator sector angle (0,) was investi- 
gated at a sector location (XJRi = 3.0), for a fixed spacing angle, ct, t = 15.0°. In the second 
grouping, i.e. Configs. 5 through 10, the influence of vortex generator sector location (XJR,) was 
investigated for a fixed spacing angle (a,*) and sector angle (0.) of 15.0° and 127.5° respectively. 
The third grouping, i.e. Configs. 1 1 through 14, investigated the effect of vortex generator spacing 
angle (a,,), keeping the generator sector angle (0,) fixed at 127.5°, and the sector location fixed 
at XJR, = 5.0. In the forth grouping, i.e. Configs. 15 through 18, the effect of generator sector 
angle was repeated, but at the axial sector location XJRi = 5i0. The last grouping of configura- 
tions was intended to show the effect of Reynold's number on the installed vortex generator 
performance. 

Distortion and Vortex Interactions with Secondary Flow 

In order to design a system of generators for vortex flow control which remains effective 
over a wide range of operating conditions, it becomes important to understand the effect of 
Reynold's number on engine face distortion as measured by standard distortion descriptors. 
Presented in Figs. 14 and 15 is the effect of Reynold's number on the radial pressure and 
60°-sector circumferential pressure ring distortion respectively, as defined by Equs. (4.0) and (5.0). 
The advantage of expressing distortion in terms of a ring descriptor is that the distribution of 
both radial and circumferential distortion can be quantified and tracked as a function of engine 
face radius. It is evident that the low energy region associated with the induced vortex pair en- 
larges between a Reynolds number of 16. Ox 10* and 4.0x10* (Fig. 14), and that this enlargement 
was accompanied by an increase in the peak 60°-sector circumferential ring distortion from 0.045 
to 0.082 (Fig. 15). In addition, peak 60°-sector circumferential distortion is moving radial to- 
wards the engine face centerline as the Reynold's number decreases. 

Presented in Figs. 16 through 20 are various aspects of vortex interactions with the sec- 
ondary flow generated within the 727/JT8D-100 center inlet duct induced by the vortex generator 
installations defined in Tables I through III. All of these interaction are described in terms of the 
60°-sector circumferential ring distortion descriptor, which is a function of the engine face radius. 
The effect of sector angle on engine face circumferential distortion (Configs. 1 thru 4) is presented 
in Fig. 16, while the influence of sector location (Configs. 5 thru 10) is presented in Fig. 17. In 
each of these cases, the vortex generator height and generator vane angle of attack were held 
fixed, along with the number of vortex generators and spacing angle. Thus, for this installation 
of vortex generators, there is an optimum axial sector location, between 5.0 and 6.0,. which will 
minimizes the 60°-sector circumferential distortion. This phenomenon will be discussed in greater 
detail in the next section. 

Choosing the axial sector location of 5.0 as the optimum station for this installation of 
co-rotating vortex generators, Fig. IS and 19 present the effects of spacing angle and sector angle 
on the 60°-sector circumferential ring distortion as a function of engine face radius. These in- 
stallations of generators are Configs. 1 1 thru 14 and 15 thru 18. 
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The effect of Reynold's number on the installed performance of vortex generator Config. 
8, as measured by the 60°-sector circumferential pressure ring distortion descriptor, is presented 
in Fig. 20. The vortex generators installation was designed at a Reynold's number of 1 2 -OjtIO 6 , 
and it is quite apparent that the minimum distortion occurs at this condition. Note also the sig- 
nificant difference in both the level and shape of the the distortion signature as a function of en- 
gine face radius. The Reynolds number effect on installed vortex generator performance will be 
examined in greater detail in the following section. 

Design Characteristics of Co-Rotating Generators for Vortex Flow Control 

Figs. 22 through 33 present a summary of the installed performance of vortex generator 
configurations 1 through 22 in terms of the engine face average total pressure recovery ( PtJPh , ) 
and the "peak" or maximum 60°-sector circumferential ring distortion descriptor that were pre- 
sented in Figs. 14 thru 20. Also included on these figures are the engine face total pressure re- 
covery' maps for each of the cases considered in this study. The engine face total pressure 
recovery (Pt.JPto) is presented both as a mass flow averaged value over the cross-sectional com- 
putational mesh, which was composed of 99 radial and 49 circumferential points in the half plane, 
and as area averaged value over a standard 40-probe rake. The standard 40-probe rake charac- 
teristics were included in this this study in order to determine whether this rake is sufficant to be 
used for vortex generator experiments. Bear in mind that because the computations were per- 
formed in the half-plane, the 40-probe standard rake compares with 9,702 computational mesh 
points in the plane of the engine face. A schmatic diagram of a 40-probe standard rake showing 
the individual probe locations is presented in Fig. 21. 

The effect of Reynolds number on engine face total pressure recovery and peak 60°-sector 
circumferential pressure ring distortion is presented in Figs. 22 and 23 for the baseline config- 
uration without vortex flow control. There is a significant decrease in the total pressure recovery', 
from 0.982 to 0.979, and an increase in the maximum circumferential pressure ring distortion, 
from 0.045 to 0.087, over the Reynolds number range from 16.0 jc 10® to 4.0x10*. Note also that # 
area averaged engine face total pressure recovery using a standard 40-probe rake poorly repres- 
ents both the level and trends of this Reynolds number phenomenon. 

Presented in Figs. 24 and 25 is the effects of vortex generator sector angle on engine face 
total pressure recovery and maximum 60°-sector circumferential pressure ring distortion at a axial 
sector location of 3.0, and generator spacing angle of 15.0°. As the number of vortex generators 
is increased, at a constant spacing angle, the sector angle angle increases according to Equ. (7), 
and this has the effect of decreasing the peak 60°-sector circumferential ring distortion as shown 
in Fig. 25. Since there are a discrete number of generators that can be placed around the inside 
periphery of the inlet duct, a minimum distortion will occure at a sector angle of 180°. 

Figs. 26 and 27 presents the effect of vortex generator sector location on the average en- 
gine face total pressure recovery' and the peak 60°-sector circumferential pressure ring distortion. 
For the installation of co-rotating vortex generators defined by Configs. 5 thru 10 in Tabel I, the 
optimum axial location lies between 5.0 and 6.0, and the overall circumferential distortion within 
this range of axial positions is less than 0.01. Thus, for a given co-rotating vortex generator in- 
stallation, there exists an axial position within the inlet duct which provides a minimum engine 
face circumferential distortion. 

Presented in Figs. 28 and 29 are the engine face average total pressure recovery and the 
maximum 60°-sector circumferential ring distortion for vortex generator Configs. 1 1 through 14. 
As the vortex generators spacing angle increases, (by decreasing the number of vortex generators 
while keeping the sector angle fixed at 127.5°), the maximum 60°-sector circumferential ring dis- 
tortion increases very rapidly above 20.0°. Thus, there exists a maximum spacing of vortex gen- 
erator blades above which the engine face circumferential distortion increases very rapidly. This 
suggests a dramatically different design guideline for vortex generator installations from that in 
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Ref. 6, but bear in mind that the effectiveness parameter used by Percy was retention of the in- 
dividual vortex identities downstream of the generator blades as measured on a flat plate, while 
the effectiveness indicator used in this study was the maximum 60°-sector circumferential dis- 
tortion. Increasing the vortex generator spacing angle does indeed increase the retention of the 
individual vortex identities, as can clearly be seen from the series of compressor face total pres- 
sure maps presentd in Fig. 29, however, while this design guideline is good for supressing local 
flow separation, it actually increases the circumferential distortion index. There was a primary 
flow separation in the baseline 727/JT8D-100 inlet, both in the experimental and computational 
duct, and this separation was associated with vortex liftoff. Each of the vortex generators in- 
stallations presented in Fig. 29 supressed this local flow separation. However, if a circumferential 
distortion descriptor is used as the measure of effectiveness, then it is desirable to have closer 
spacing between vortex generators than would normally be used to supress local flow separation. 

The effect of the vortex generator sector angle on the maximum 60°-sect or circumferential 
pressure ring distortion at an axial sector location of 5.0, (Configs. 15 thru 18 in Table II), is 
presented in Figs. 30 and 31. The relationship between number of vortex generators and the 
sector angle for a fixed spacing angle is again given by Equ. (6.0). In comparison with a similar 
study done at the axial sector location of 3.0 Fig. 25, it is clear that this installation of co-rotating 
generators is best located at an axial location of 5.0. Had another installation of generators with 
different physical dimensions been chosen, then a different optimum axial sector location would 
have resulted. 

The effect of Reynolds number on engine face recovery and distortion for the vortex gen- 
erator Config. 8 is presented in Figs. 32 and 33. For this installation of vortex generators, the 
average total pressure recovery level (Pi.r/Pk) and maximum 60°-sector circumferential pressure 
ring distortion index remain reasonably level between the Reynolds numbers of 1 6.0jc 1 0 6 and 
8-OjrlO 6 . For Reynolds number less than 8-O.xlO 6 the flow at the engine face "breaks" down and 
the distortion increases very rapidly. The systematic and continuous nature of the flow field 
breakdown can be seen in the engine face total pressure recovery maps presented in Fig. 33. In- 
stalled vortex generator performance, as measured by standard engine face circumferential dis- 
tortion descriptors, is sensitive to Reynolds number and thereby the generator scale, i.e. the ratio 
of generator blade height to local boundary layer thickness. Installations of co-rotating vortex 
generators work well in terms of minimizing engine face distortion within a limited range of gen- 
erator scales. Clearly this effect is not well understood since it may be influenced by the structure 
and strength of the induced secondary flow. Consequently, this is an area for future research. 
However, it is important to recognize that the design of the vortex generator installation is a 
point design, and all other conditions are off-design. 

Comparison with the 727/JT8D-100 Comformation Test Data 

The Phase II low speed conformation tests on the 727/JT8D-100 center inlet duct per- 
formed by Kaldschmidt, Svltedo, and Ting Ref. 7, were sponsored by NASA Lewis Research 
Center and performed in the Boeing 9 ft. x 9 ft. low speed wind tunnel between July and No- 
vember of 1973. This series of experiments were performed as a conformation test for a re- 
designed center inlet duct for the 727-100 series aircraft, and they w^ere never meant to be used 
as validation for computationsl codes. As such, the information presented in Ref. 7 is incomplete 
for code validation, however the data is valuable enough to attempt a comparison with compu- 
tations using RNS3D. It must be understood that much of the input information was estimated 
or gleaned from figures or photographs presented in Ref. 7, and do not represent precise values. 

The experimental and computed engine face total pressure contours for the 727 JTSD-100 
baseline duct without vortex generators is shown in Fig. 34. The experimental model had an 
engine nose dome of elliptical shape of aspect ratio 2.0, while the calculations were performed 
without the dome for simplicity. In general, the comparision between analysis and experiment 
is quite good, bearing in mind that the exact tunnel conditions were not specified in Ref. 7, nor 
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is the influence of the nose dome understood. The effect of the nose dome on the engine face 
total pressure recovery map and distortion level has never been systematically studied, but some 
performance advantage might be gained since the circumferential distortion descriptors are very 
sensitive to both level and distribution. 

Fig. 35 present a comparision between the engine face total pressure map with the "best" 
experimental vortex generator installation (Config. 12 in Ref. 7), and RNS3D results with the 
generator pattern which can be gleaned to be most similar to the experimental installation 
(Config. 3 in Table I). Both vortex generator installations were composed of seven pairs of co- 
rotating generators assumed to be symmetrically located about the 180° azimuthal reference, and 
positioned at an axial location of 3.0. However, the experimental generator installation also 
contained three pairs of counter-rotating generators located symmetrically about the 0.0° 
azimuthal reference, and positioned just upstream of the engine face, while the calculations did 
not include any counter-rotating generators. In general, calculations of the vortex generator flow 
field, and in particular the engine face total pressure distortion map, shows remarkable agreement 
with the measured engine face recovery map. 

Presented in Fig. 36 is a comparison between the engine face recovery maps resulting from 
the experimental vortex generator installation judged the most effective (Config. 12 in Ref. 7) and 
the computational generator installation judged most effective (Config. 15 in Table II). The 
computational vortex generator installation was composed of eleven pairs of co-rotating genera- 
tors located at an axial sector location of 5.0. While the computational vortex generator instal- 
lation had one additional pair of generators over the experimental configuration, the additional 
loss associated with this pair will be remarkable small as indicated in Table IV. The measured 
loss level resulting from ten pairs of vortex generators is 0.002 as presented in Ref. 7, thus the loss 
level associated with the additional pair will be very small. The loss levels associated with prop- 
erty design vortex generators installations are in general very small, and thus they represent a very 
good method to manage engine face distortion by restructuring the development of the secondary 
flow that is induced within three-dimensional curved inlet ducts. 

CONCLUSIONS 

The present results provide a validation of the initial value space-marching 3D RNS pro- 
cedure and demonstrate accurate predictions of the compressor face flow field, with a separation 
present in the inlet duct as well as when vortex generators were installed to suppress separation. 
The computing time on the CRAY XMP (i.e., CPU = 6.5 min. for 2.28 x 10 s grid points) for 
University of Tennessee diffusing S-duct configuration for both the baseline and the the case with 
three pairs of counter-rotating vortex generators, is sufficiently rapid for routine use in an anal- 
ysis and design environment. 

This paper demonstrates that the three-dimensional Reduced Navier-Stokes code RNS3D 
can be used very effectively to develop a vortex generator installation for the purpose of mini- 
mizing engine face distortion by restructuring the development of secondary flow induced in 
typical three-dimensional curved inlet ducts. For a given installation of vortex generators, there 
exists an axial location within the inlet duct which will minimize engine face distortion, and within 
this installation, there exist a maximum spacing of generators above which the distortion in- 
creases rapidly. The installed performance of vortex generators is sensitive to Reynolds number 
or generator scale, i.e. the ratio of blade height to local or average boundary layer thickness. 
Consequently, there is a limited range of generator scales that serve to minimize engine face dis- 
tortion. Therefore, the design of vortex generator installation is a point design and all other 
conditions are off-design. In general, the loss level associated with a properly designed installa- 
tion of vortex generators is very small, and thus they represent a very good method to manage 
engine face distortion. This study also shows that the vortex strength, generator scale, and sec- 
ondary flow field structure have a complicated and interrelated influence on the engine face dis- 
tortion, over and above the influence of the initial arrangement of generators. 
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Figure (2) - Computation grid for the Univ. Tennessee diffusing S-duct. 
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Figure (5) - Reduced Navier Stokes (RNS) solution showing secondary 
flow structure from vortex generator model, X/Z>,-= 0.18. 



Figure (6) - Reduced Navier Stokes (RNS) solution showing limiting 
streamline signature of vortex generator region. 


17 








Figure (7) - Reduced Navier Stokes (RNS) solution showing limiting 

streamline signature of separation region within Univ. Tennessee 
diffusing S-duct. 



Figure (8) - Surface oil floM patterns showing separation region 
within Univ. Tennessee diffusing S-duct. 
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Figure (9) - Geometry definition for the 727/JT8D-100 center inlet. 



Figure (10) - Computation grid for the 727/JT8D-100 center inlet. 
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Figure (1 1) - Geometry definition of co-rotating vortex generators. 



Figure (12) - Nomenclature used for vortex generator positioning. 
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Xyg/Ri = 4.0 


Figure (13) - Axial locations of the vortex generator sector regions. 


Config. No. X,gJRi ii, x ac,. x 0, Reynolds No. 

1 3.0 

11 15.0 157.5 12.0x10* 

2 3.0 

9 15.0 127.5 12.0x10* 

3 3.0 

7 15.0 97.5 I2.0.vl0* 

4 .3.0 

5 15.0 67.5 12.0x10* 

5 2.0 

9 15.0 127.5 12.0x10* ) 

6 3.0 

9 15.0 127.5 12.0x10* | 

7 4.0 

9 15.0 127.5 12.0x10* ; 

8 5.0 

9 15.0 127.5 12.0x10* l 

9 6.0 

9 15.0 127.5 12.0x10* % 

10 7.0 

9 15.0 127.5 12.0x10* \ 


Table 1 - 727/JT8D-100 center inlet duct vortex generator configurations. 
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Table II - 727/JT8D-100 center inlet duct vortex generator configurations. 



Table III -727/JT8D-100 center inlet duct vortex generator configurations. 
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Figure (14) - Effect of Reynolds number on radial pressure ring 
distortion without vortex flow control. 



Figure (15) - Effect of Reynolds number on the 60°-sector circumferential 
pressure ring distortion without vortex flow control. 
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Figure (16) - 


Effect of vortex control sector angle (0,) on the 60°- sector 
circunferential pressure ring distortion at XJRi = 3.0. 



Figure (17) - Effect of vortex control sector location (XJR) on the 60°-sector 
* circumferential pressure ring distortion. 







arc M f min 



Dimensionless Engine Face Radius, S ef 


Figure (18) - Effect of vortex control spacing angle (a„) on the 60°-sector 
circumferential ring distortion at X r JR,- = 5.0. 



Dimensionless Engine Face Radius, 


Figure (19) - Effect of vortex control sector angle (0,) on the 60°-sector 
circumferential pressure ring distortion at X UR, = 5.0. 



Re = S.O.vlO 6 
Re = 6.0.vl0 f ' 
/?c = 8-O.vlO 6 
Re = 12.0.vl0 6 
-/?<? = 16.0.vl0 6 


IB 


Dimensionless Engine Fnee Radius, S t f 

Figure (20) - Effect of Reynolds number on the ^-sector circumferential 
pressure ring distortion at A, ; / a, 5.u. 



Figure (21) - Standard 40-probe rake geometry and probe locations. 
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Reynolds Number, Re x 10 -< 


Figure (22) - Effect of Reynolds number on the engine face total 

pressure recover) (PtJPt 0 ) without vortex flow control. 



Reynolds Number, Re x 10~‘ 


Figure (23) - Effect of Reynolds number on the maximum 60°-sector 

circumferential pressure ring distortion without vortex flow control. 
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Total Pressure Recovery, PtjPt 0 



Figure (24) 


Effect of vortex generatojjiector angle (0.) on the engine face 
total pressure recovery (Pt^Pto) at X,,ID, = 3.0. 



Vortex Generator Sector Angle, 0 , , degrees 

Figure (25) - Effect of vortex generator sector angle (^ on the maximum 
8 1 60 = -sector circumferential pressure ring distortion at A ,,/a, 3. . 
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Total Pressure Recovery, N t/ fPt n 



Figure (26) - Effect of vortex generator sector location ( X,JR ,) on the engine 
face total pressure recovery' (Pt e/ ]Pt 0 ). 



Axial Distance, X^/R, 


Figure (27) - Effect of vortex generator sector location (X.JR) on the maximum 
60°-sector circumferential pressure ring distortion. 
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Vortex Generator Spacing Angle, , degrees 


Figure (28) - Effect of vortex generatorspacing angle (a,,) on the engine face 
total pressure recovery ( PtjPt 0 ). 



Vortex Generator Spacing Angle, a,, , degrees 


Figure (29) - Effect of vortex generator spacing angle (a,,) on the maximum 
60°-sector circumferential pressure ring distortion. 
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EXPERIMENTAL DATA ANALYSIS 

Figure (34) - Engine face total pressure contours for the 727/JT8D-100 
inlet duct without vortex generators. 


3 Pairs Counter-Rotating Generators No Counter-Rotating Generators 



7 Pairs Co-Rotating Generators 7 Pairs Co-Rotating Generators 


EXPERIMENTAL DATA, CONFIG. 12 ANALYSIS, CONFIG. 3 

Figure (35) - Engine face total pressure contours for the 727/JT8D-I00 
inlet duct with vortex generators. 
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Figure (30) - 


Effect of vortex generatoj^sector angle (0,) on the engine face 
total pressure recover)' 0 ) at X,,ID, = 5.0. 



Figure (31) - 


Effect of vortex generator sector angle (8.) on the 

60°-sector circumferential pressure ring distortion at A,,/*, - 5.u 
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c Total Pressure I 



•e (32) - Effect of Reynolds number on the engine face total 
pressure recovery’ (Pt ef fPt 0 ) at XJRi = 5.0 



Reynolds Number, Re x 10 - * 


Figure (33) - Effect of Reynolds number on the maximum 60°-sector 
circumferential pressure ring distortion at A = 5.0 



3 Pairs Counter-Rotating Generators No Counter-Rotating Generators 



7 Pairs Co-Rotating Generators 11 Pairs Co-Rotating Generators 

EXPERIMENTAL DATA, CONFIG. 12 ANALYSIS, CONFIG. 15 

Figure (36) - Engine face total pressure contours for the 727/JT8D-100 
b inlet duct with vortex generators. 



Table IV- Summary of experimental and computed engine face total 

pressure recovery ( PtJPto ) for the 727/JT8D-100 inlet duct. 
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